Dephasing of Josephson oscillations between two coupled Bose-Einstein condensates 



O 

oo 



Pierre Villain 1 ' 2 and Maciej Lewenstein 2 
1 Commissariat à L'Energie Atomique, DSM/DRECAM/SPAM 
Centre d'Etudes de Saclay, 91191 Gif-sur-Yvette, France 
2 Institut für Theoretische Physik, Universitat Hannover 
Appelstr.2, 30167 Hannover, Deutschland 

We study the dynamics of the relative phase between two atòmic Bose-Eistein condensates 
coupled via col·lisions and via a Josephson-like coupling. We derive the equations of the motion 
for the relative phase and the relative number operators from the second quantized Hamiltonian of 
the system using a quantum fïeld theoretical approach. We distinguish the cases in which the two 
condensates are in the same trap or in two different traps and study the influence of this difference 
on the first order correlation function of atòmic fields. In identical traps this function does not 
on ■ undergo dephasing. We calculate the dephasing time for the case of different traps. 
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^ 1 I. INTRODUCTION 

Since the observation of the Bose-Einstein condensation (BEC) in trapped cold alkali gases a lot of attention 

7-H ' has been devoted to explore the coherence properties of the condensates. In particular, in a recent experiment, the 
interference pattern between the two overlapping condensates has been measured Q, in very good agreement with 
subsequent theoretical calculations In the context of interference, the static and dynamic properties of the phase 
of the condensate are of major importance, and for this reason they have been the object of many theoretical studies in 
the recent years. First of all, the question of very existence of the phase has been put forward. The Standard theory of 
BEC that is Bogoliubov-Hartree-Fock theory Q assumes implicitely that the condensation in the weakly interacting 
Bose gas occurs into a coherent state of the lowest energy mode of the system. Such a state admits a defined phase 
ON 1 and amplitude, but exhibits inevitable fluctuations of both these quantities. In particular, it assumes fluctuations of 
■ the number of àtoms in the condensate. On the other hand, if the atom number in a system is conserved, then it is 
Oh! clear that at very low temperature the condensate will be in a Fock, or in another words number state, which excludes 
1 any atom number fluctuations 0. Moreover, it is not possible to assign a definite phase to such Fock state. 

Several authors have shown, however, that two overlapping condensates show similar interference patterns, inde- 
pendently of whether they arc in the coherent or Fock states P-[ïo|. The reason is that repetitive measurements on 
two condensates lead to a build-up of the relative atom number fluctuations, and consequently, to a determination of 
L* ' a phase. These results, stimulated by older discussions on the existence of the phase of the condensate jn]], has led 
. £h [ to a conclusion that many properties of systems with large number of partides can be described either way, i.e. using 
^ ■ coherent states and their phases, or avoiding it [Q . 

The qüestions what is the character of the quantum state into which the condensation occurs, and what are its 
properties are nevertheless of major importance. In particular, recently several authors have pointed out that the 
relative phase of two condensates, if established, must necessarily undergo collapse and revivals due to non-linear 
dynamics of interactions between àtoms ]Ï3|-p^| . The dynamics of the phase refiects the quantum properties of the 
state of the system, and in particular the collapse time is inversely proportional to the dispersion of the relative 
number of àtoms in the two condensates system. 

Originally, the calculation of the phase properties has been done with a single condensate In particular, in 

Ref. Q we have developed a fully quantum fïeld theoretical method to describe the phase dynamics and to calculate 
the collapse, assuming a breaking of the U(l) (phase) symmetry of the system. Such a symmetry breaking admits 
the existence of fluctuations of the phase and of the number of àtoms in the condensate. The results concerning the 
dynamics of an absolute phase and absolute atom number for a single condensate can be straightforwardly applied to 
the case of a relative phase and atom number for double, or múltiple condensates in stationary states [|Ï4|jÏ7]]. 

In fact, in Ref. [^4| we examined two aspects of the phase diffusion or collapse in the stationary state of two 
condensates, that overlap and may even interact through elàstic col·lisions. First, we have calculated the properties of 
the quantum state of such system after a preparation phase, during which the condensates were coupled via coherent 
Josephson-like coupling. Such coupling leads to fluctuations of the relative atom number, and determines a well 
defined relative phase between the condensates in the ground (stationary) state of the system. After the preparation 
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period, the coupling between the condensates is turned off, and the condensates may evolve freely, or remain coupled 
through the phase and relative number preserving col·lisions. In both cases the memory of the prepared phase is 
lost after the collapse time which is now given by a sum of two contributions of the two condensates, each inversely 
proportional to the dispersion of the corresponding atom number. 

It should be stressed that although the revival times are typically beyond the reach of present experiments, the 
collapse phenomenon and its influcnce on the phase memory could be in principle measured. So far, the JILA group 
has realized two overlapping condensates of rubidium àtoms in (F = 1, tuf — 1) and (F — 2, mp — 2) in a magnètic 
trap |Ï8| . The overlap in such a system is only partial, but nevertheless allows in principle to transfer àtoms coherently 
from one state to another using appropiate sequence of radiofrequeney and microwave transitions (alternatively, far 
resonant stimulated Raman transition could be used). In a more recent paper |Ï9| , the JILA group reported the 
results of experiments aiming in testing the robustness of the relative phase between two overlapping condensates. 
Their results show that the random phase shift due to the coupling with the environment is not significant on a during 
of the order of lOOms. Conversely the effects of the quantum phase fluctuations were negligible on this time seale. 
Namely, in the conditions of the JILA experiments (5.10 5 àtoms in a trap with axial frequeney 59 Hz) we estimate 
the diffusion time to be of the order of 0.5 s. So that the time linear shift due to quantum fluctuations could not be 
noticeable. It is to be noticed that these effects could be seen with a condensate formed of less àtoms. In the recent 
experiment of the MIT group p(J three co-existing condensates in the ground state of sodium (F = 1; = — 1,0,1) 
have been created by transferring the (F = l,m = —1) condensate from the magnètic to a purely optical dipole 
trap. In such trap àtoms in all Zeeman states remain trapped, and moreover, feel the same trapping potential. Using 
radiofrequeney transition one can transfer àtoms from one state to another. 

Both above mentioned experiments will allow thus to study coherent Josephson oscillations between the two, or 
more condensates. This problem has been already studied theoretically, and the corresponding nonlinear Schròdinger 
equation for the condensate wave functions has been investigated numerically pï| . To our knowledge the problem 
of quantum fluctuations around dynamical solutions of the coupled Gross-Pitaevskii equations (i.e. the problem of 
coupled time dependent Bogoliubov-de Gennes equations) have not been studied so far. Some attempts to study 
the quantum dynamics have only been made for the simplified two mode models p2| . In particular, the problem of 
dephasing of the Josephson oscillations has not been solved using a full quantum field theoretical approach. Thcrc 
are two physical aspeets of such dephasing: first, it oceurs due the intrinsic nonlinearity of the Josephson oscillations, 
that is the fact that oscillation frequeney depends, in general, on the inital state of the system. Since any quantum 
state exhibits nccessarily quantum fluctuations, it is to expect that the frequeney of the corresponding Josephson 
oscillations would exhibit some fluctuations too, leading thus to dephasing. Second, the quantum nature of the 
system may lead to a spreading of the distribution of the relative phase between the two condensates simultaneously 
with the Josephson oscillations. In this paper we investigate the interplay between those two effects. 

In section 2 we describe our model of two interacting condensates, and derive the Heisenberg equations of the 
atòmic field operators. Then, using a Bogoliubov approach W, we obtain a set of coupled Gross-Pitaevskii equations 
for the mean vàlues of the fields operators, and the coupled time-dependent Bogoliubov-de Gennes equations for the 
quantum fluctuating parts of the field operators. In section 3 we introduce the operators describing the fluctuations of 
the total and relative number of àtoms and those describing the total and relative phase, respectively. In section 4 the 
dynamical equations for these operators are derived using a rotating wave approximation, which allows to decouple 
the slow evolutions of these operators from the relatively faster evolution of the quasi-particle operators. In section 
5 we solve analytically the dynamics of the phase and atom number operators in the case when the two condensates 
are trapped by the same potential (symmctric case). We show that the first order correlation function of the atòmic 
fields of the two condensates does not present any dephasing in time. In section 6 we study the case of two different 
trap potentials (the assymetric case), and solve the equations perturbatively with respect to the relative potential 
difference. Our study shows that there exists a dephasing, characterized by a time seale whose dependance with the 
different parameters of the problem is calculated. We conclude in section 7. 

II. DESCRIPTION OF THE MODEL 

Our system is composed of two condensates consisting of àtoms in two different internal states A and B in two 
overlapping harmònic traps characterized by the different freqüències u>a and lob- The difference exists in general in a 
magnètic trap since different Zeeman states react differently to the magnètic field; nevertheless, the freqüències might 
be the same if the trap is purely optical such as recently used by the MIT group p0[ . At very low temperature the 
interaction between the àtoms is well modeled by a contact interaction potential uo8(r — r ') with uq = 4:nh 2 asc/ m 
where ase is the scattering lenght of the collision process. In our case we consider three different scattering processes: 
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the A — A, B — B and A — B col·lisions, and thus three different scattering lenghts. For the case of the Rubidium 
as in the JILA experiment we can neglect the differences between these scattering lenghts and assume that the 
interaction potential is the same for every atom regardless of its internal state. In the following we will therefore write 
uaa = ubb = u ab = uq. We assume also a coherent Josephson like coupling between the A and B àtoms which can 
be realized using a làser (Raman) transition between the two states A and B. In order to make this process efficient 
it is necessary to have a non-negligible overlap between the two condensates. This last condition can be achieved in 
the JILA case |ïï| (which will be our reference experiment for the general study of the two differents traps), as well 
as in the optical trap experiment of MIT . 

The second quantizcd Hamiltonian of the system is: 



H 



dr^ A (f, t) 



.ií_V 2 + V^(f) 
2M Ay ' 



«0 



+ y / dr$>\{r, m\{r, t)V A (f, t)V A (r, t) 



df^ B (f,t) 



drÍ> B (r, t)¥ B {f, t)Ü B (? , t)Ü B (r, i) 



+ u I df¥ A (f, t)¥ B {f, í)* B (r , t)V A {r , t) 
-h\ ídr \$ B (f, t)Ü A {f, t) + ¥ A {f, t)Ü B (r, i) , 
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where ^> A ,B(r,t) denote the atòmic field operators fullfilling the Standard (bosonic) commutation relations: 

^ A , B (f,t),¥ AiB (r',t)\ =5(?-f'). (2) 

The last term in the Hamiltonian (|ï|) represents the Josephson like coupling which allows exchange of àtoms from one 
state to another. The frequency A characterizes the strenght of the coupling between the two atòmic states and its 
value can be experimentally controllcd. 

From the Hamiltonian (fil) we derive the Heisenberg equations for the atòmic field operators: 



= [Ú A (r,t),H] 



$ A + U Q Ú A Ú A ^A + U Í/ B ^B^A - h\u Í/ B , 



and 



ih 



dt 



[* B (f,t),H] 



2M 



V B (r) 



(3) 



(4) 



We use then the well-known time dependent Bogoliubov method to find the dynamics of the system at T — 0. To 
this aim we assume that the atòmic fields operators can be written as : 

(5) 
(6) 

where N denotes the total number of àtoms, \/Nij) A {f, t) and y/Ni/j B (r,t) are in fact the wave functions of the A 
and B condensates whereas the S^a,b denote the fluctuating parts of the operators, respectively. By setting these 
expressions into the Eqns (||) and (Q) and identifying the parts of zero and first order respect to the fluctuations, we 
find the following set of coupled equations: 



V A (r,t) = VNtp A {r,t) +Stf) A (f,t), 
* B (f , t) = ^ÍN^ B {r, t) + 5ip B (r, t), 
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ih = C A ij) A (r,t) +u (p A + p B )ijj A (f,t) - h\ip B (r,t), 



iti 



dt 

dip B (r,t) 
dt 



C B 4> B (r,t) +u (pa + p B )tJj B (r,t) - h\ip A (r,t), 



(7a) 
(7b) 



for the zero order part. Here C A:B = ~jm^ 2 + Va,b(t) , while p A .s — N\t(j AtB (f , t)\ 2 are the densities of the A 
and B àtoms. In the following we will denote by p — p A + ps the total density of àtoms. Moreover we assume the 
following normalisation conditions: Jd 3 r N\tp AtB (r ,t)\ 2 = N AtB (t) where N A and Nb denote the number of the A 
and B àtoms. Obviously in this case we have J d 3 f [\ip A (r, t)\ 2 + \ip B (r, t)\ 2 ] — 1 so that the total number of àtoms 
is fixed to N. Eqs. ([7a 7b) constitute a system of coupled Gross-Pitaevskii equations. 
In the first order in the fiuctuating parts we find: 



ih = [C A + u p\ òip A + Nu ip A 



ih 



dt 

dStp B (f,t) 
dt 



ip A 6ip\ + ip* A 5ip A 



= [C B + u p] 54> B + Nuqi^b tPbSvb + ^*B^ 



Nu ip A 

- Nu 1pB 



^bH b + ^B^B 



hXSi/jB, 



(8b) 



These equations in turn constitute a system of coupled Bogoliubov-de Gennes equations ( |6|,|23|]). As the Hamiltonian 
is invariant under the permutation of the A and B àtoms, it is convenient to introduce the following atòmic field 
operators: 



Ú + = e - iM [Ü A + Ú B ]/V2, 



= e lM fi> A - *i 

Using the previous evolution equations we can easily derive the equations for the new wave functions tp+ and ip 

.dip+{r,t) 



(9a) 
(9b) 



ih 



dt 
dt 



= [C + u p]ip + + SViJj-, 
= [C + u oP ] ^_ + ÍV*V+> 



(10a) 
(10b) 



where C 



2M 



V 2 + V m {r) , SV = SVe- 2lXt , V m = (V A + V B )/2, SV = (V A - V B )/2. Notice that the unitary 



character of the transformation (|9a^b) leaves the total density of àtoms invariant: p — p A + pb = P+ + P-, with 
p± = N\ip±\ 2 . For the fiuctuating parts of the operators we obtain : 



ih 



dt 

dSijj-(r, t) 
~dt 



[L + u (p+ + 2/j_)] 8ip- + SV*StP + + Nu Tp 2 _5ipl + Nuq^-^+Si/jI + ^1Si/j + ). 



(lla) 
(11b) 



III. TOTAL ATOM NUMBER AND PHASE 



In this section we introduce the operators describing the total atom number and phase fluctuations. We are 
following the method previously developed in refs. JÏ3|,[Ï4[ to study the same system, but in the static case. We begin 
by introducing the following hermitian operators: 

P+= J tí>%6$+ + iP+ó^ + )d 3 r, (12a) 

P_ = J \ilj*_5i)-+^-5^_)d 3 r. (12b) 

These operators describe, r espe ctively, the fluctuations of the number of àtoms in the new states ip+ and Using 

llh| ) we obtain the evolution equations of this two operators in a straighforward way: 



the equations (10a,10bj,lla, 
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ih 



dP + 
~~dT 



SV 
dP. 



d 3 r- 



5V* 



d 3 r, 



-ih- 



dt dt 
We then construct the operators: 



+ =ih (^+ - «M^)^, 



ih 



f_5i{>- - (j)-6ij)l)d 3 r, 



for which thc commutation rclations read: 



Q+,P+ = 2ifíReJ + , 



where 



Q-,P- 

J+ = 
J_ = 



= 2ihRe J_ 

cf)* + ip + d 3 r, 
(j)*_ïP-d 3 r. 



(13a) 
(13b) 

(14a) 
(14b) 

(15) 
(16) 

(17a) 
(17b) 



In the following wc will dcnote Re J± by j± . 

The operators (14a, 14b ) up to a normalisation factor are conjugated to the operators (12a, 12bj ) and describe the 
fluctuations of the phases of the new atòmic field operators. In order to study their dynamics, we have to find 



the evolution equations for the functions 4>+ an d 



For this purpose we introduce the operators describing the 



fluctuations of the total number of àtoms and the total phase of the system: 

Ptot = (P++P-)/V2, 



(Q+ + Q-)/V2 



The latter operator fulfills ( 110): 



Qtot-, Pu 



and 



dQ t 



dt 



= NÜPt, 



(18) 
(19) 



(20) 



(21) 



where the coefficient ü will be define in the following. The commutation relation (J2Q) imposes the following rclation 



Re J+ + Re J_ = 1 

and the operator evolution equation (|2^) gives the evolution equations for (f> + and </>_ : 

.d</>+(r,t) 



(22) 

(23a) 
(23b) 

In order to define in a unique way the functions 0+ and </>_ and the coefficient ü (indeed the functions <j>'± = 4>± + iatp± 
where a is a real coefficient are solutions of the above equations too) we impose a more restrictive normalisation 
condition by requiring additionally: 



ih- 
ih- 



dt 
dt 



= [£ + u a (2p + + p_)] (j)+ + Nuotp^l + Nu ip+(ip-^*_ +ipl<f>-)- Nüip + + 5V(j)-, 
= [£ + u ( P+ + 2pJ)\ ct>- + Nu ij 2 _<p*_ + Nu ip~{ÍJ + (f>* + + ip* + (/)+) - Nüip- + 5V*<j>+. 



ò* + tP+ + (t)^*_)d 3 r = 1, 



(24) 
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id est 



J+ + J_ = 1. 

It is interesting to note that there exists a direct relation between <j)± and ip±, namely: 

0± = 2tffe-Í/(t)V>±], 



dN h J 



with 



and 



9^+ 3 



<9A 



<97V 



ü = -2^L = ~2h— 
dt dt 



Im / [V> 



difj 



+ 



+ <9a 1 r ~ sa 

The exact Heisenberg equations for the phase operators have the form : 



dQ 4 
dt 



dQ- 
dt 



NÜP+ + Nu / (ip+(f>l + ipl(f> + )(ip-Si)L + iP*_5ip^)d 3 r - Nu / (ip-(f>*_ + ip*_<f>-)(ip + Si>l + ^*+H+)d 3 r 



(25) 



(26) 



(27) 



(28) 



+ / 6V(<pl5ip- - cj)-Sipl)d 3 r + I 8V*(^+Si/?l - <t>*_54,+ )d ó r, 



(29a) 



NÜP- - Nuo / (V>+</>+ + ipl(f> + )(ip-6ipl_ + ip*_6i(j-)d 3 r + Nu a / (V>-<£* + i>_ <£_)(V>+<fyí>+ + ^+<ty>+)d 3 r 



5-0- - c/>-6ip+)d 3 r - I 6V*((j) + 6ip J _ - cj)^5^ + )d 3 r. 



(29b) 



IV. EVOLUTION OF THE OPERATORS IN THE ROTATING WAVE APPROXIMATION 

It is possible to simplify the previous evolution equations by making a Rotating Wave Approximation (RWA). 
Namely, we make use the fact that the quasi-particles operators describing the elementary excitations evolve with a 
characteristic time of the order of magnitude of the period of the harmònic trap. On the othcr hand the ev olution o f 
the operators P+, P_, Q+, Q- is slower, and it is thus possible to project the evolution equations (13a,13bmd 29a , g9t ) 
on P + , P_, Q + , Q_ . Introducing the following notation: 



((j)\iP) = / 5V<j)*ip d 3 r, 



we then find in the RWA : 



dR 



+ 



dt 

dP- 
dt 

dQ+ 
dt 



dt 



-Im(^+|V-)P + 



— Im(^+|^-)P- . 
rvj- h 



\ 7+ 7- 



dP + 
~dT 



Nü- 



2Nu I 
7+ 



-Re(0 + |0_) 



P + 



[27Vu o X + Re(0+|0_)]P_ 



— Im(^ + |^_)Q H 
«7+ 



1 



Irm>+|V-)Q- 



— [2Nu l + Re(0+|0_)] P + 
7+ 



2iVu I 1 . , ' 
Nü Re(^ + 0_) 

7- 7- 



P_ 



-Llm(V>+|0-)Q+ - -^Im(0+|^)Q- 
7i7+ «7- 



(30) 

(31a) 
(31b) 

(31c) 
(31d) 
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where X = J Re(í/;-0I)Re(V'+0+)d 3 r. At this stage it is interesting to introduce the operators describing the fluctu- 
ations of the relative number of àtoms between the states and ip- and of the relative phase between these two 
states. We use the approach previously developped in 14 to dehne these operators: 



a _ÍQ ± _Q_" 

V 7+ 7- , 



2. 



(32) 
(33) 



We can then write the evolution equations of the new set of operators Ptot , Preí , Qtot, Qrel in terms of those operators: 
dPtot 



dt 

dQtot 
dt 

dPrel 

dt 

dQrel 

dt 



= 0, 



NüP to t 
1 



7+ 



Jm((f> + \ip-) 



7- 



-Im(-0+|^_) 



Prel + 



1 



Re{il> + \ip-)Q re i, 



1 



7+7- 



Nü- 



2Nu I 
7+7- 



1 



7+7- 



-Re(<^ 



7- 



— Im(0+|V>-) 
7+ 



5 rei • 



(34a) 
(34b) 
(34c) 
(34d) 



These equations constitute one of the bàsic results of the present paper. We can remark the complete decoupling 
between the evolution of the observables of the total dynamics and these of the relative dynamics. We shall now 
discuss the solutions in various limiting cases. 



V. THE SYMMETRIC CASE 



We call the symetric case the case where the two condensates A and B are in the same trap (mathematically that 
means that SV = 0). Experimentaly such situation would occur when the àtoms are confined in a purely optical trap. 
This kind of trap for Bose condensates has been recently realized by the MIT group |20| , and is very important from 
the theoritical point of view since the problem in this case can be exactly solved. In the, so-called, Thomas-Fermi 
limit one can show that the solutions of the Gross-Pitaevskii equations have the form |Ï4|,^4],^5| : 



tff (r,t) 



N 

fWZ 

N 



y/p(r,t) 



,i9_(r,t) 



where the density p is the Thomas-Fermi solution of the single condensate Gross-Pitaevskii equation 



15N 



1 



" 87rr 3 (í) V r 2 (í), 

and tq is the size of both condensates in the Thomas-Fermi approximation, fullfilling the equation: 



f + uj t (t)r 



0. 



(35a) 
(35b) 

(36) 

(37) 



The phase 8± has the form (A(t)r 2 + B±(t))/U, where the coefhcients A{t) and B +) _(í) follow the evolution equations 
previously found in Ref. E3]: 



.4 



-mr /2r , 



B± = 15iVu /87rrg, 

and N + and 7V_ are constant. We find thcn the following form for the functions 0+ and 



(38) 
(39) 
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,(0) _ o(t) 
,(o) _ a(í) 



(5P ( 4 °) 



(41) 

where the function a(í) is equal to (^f-^o^)) -1 because of the normalisation conditions ( p2| , p5| ), which requires more 

7!°"' = 7_ = i. Moreover, we have AuQOt{t) = because the functions <f>jP are solutions of the same equations 

as ip±\see H). Taking SV = in the equations ( 34a , 341 , 34c , pld ) and replacing the solutions obtained of the 
hydrodynamic limit in X we get immediatly: 

0, (42) 

NÜ^Ptot, (43) 

0, (44) 

0. (45) 

(46) 



dt 

dQtot 
dt 

dPrel 

dt 

dQrel 

dt 



The fact that Ptot and P re i are constants of the movement is related to the conservation of the total number of àtoms, 
and to the conservation of the relative number of àtoms (which is evident since both N + and iV_ are constant). 
Moreover, we conclude that the relative phase between the two states + and — is fixed (up to fluctuations constant in 
time) and does not diffuse. Only the total phase operator Qtot exhibits a linear growth in time (see Eq. ([43[)) which 
leads to a gaussian decay of the correlation function of the type < ^ A (r, t + r)^A(r, t) > as discussed in |13| . Using 
the inverse transformation : 

* A = [e íAí *+ + e- íAt *_]/\/2, (47a) 
* B = [e !Af * + - e _iAí *_]/V^, (47b) 

it is then possible to calculate the one time correlation function < ^^(f, í)$s(r, í) >. This function oscillates with 
the frcquency 2A, which means that the fluctuations of the relative phase between these two states don't show any 
effect of phase diffusion. This is a very interesting result since it immediatly implies that in weakly asymetric case, 
the phase diffusion rate will have to scale with SV. 

It is worth noticing that in this particularly simple case it is not necesary to make the transformation ( pa| , pb| ) . The 
Gross-Pitaevskii equations for ipA and i\)b can be exactly solved in the hydrodynamic limit using a selfsimilar solution 
of the form : 



( 15*k(*n * ( 1 r 2 Y AA{tV+BB{t))/ n 
\8nNrUt)J \ r 2 (í) 



tt B (F,t) = ; ,7, U--Ï7T e HAWr-+n BW ;/ B (48) 



Introducing the notations A<ï> for the relative phase $a — &b = (-Ba(í) — -Bb(í))/?í and AN for the difference of 
populations Na(í) — Ns(t), and using the equations of the motion we find: 

A *=^> (49a) 
dC 

A$ = , (49b) 

where C is given by: 



C = 2h\y/N 2 - AN 2 cos(A$). (50) 
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The Eqs. (49a, 49b| ) constitute a set of classical Hamilton equations for which the function C is a Hamiltonian, i.e. a 
constant of the motion. By straightforward integration we obtain: 

AN = Acos{2\t + <$> N() ), (51) 

where 



A = \ N 2 



4h 2 \ 2 ' 



iV 2 sin 2 (A$) + A7V 2 cos 2 (A$). (52) 

and is a constant depending of the initial conditions, namely we have $no = arccos(AiV(0)/^4). The difference 
of populations between the two condensates oscillates at the previsible frequency 2A. In the same way we find that 

AHXA 

tan(A$ - A$(0)) = — - — sin(2Aí). (53) 



VI. THE ASYMMETRIC CASE 



We want now evaluate the influence of the difference between the two potential traps on the dynamics of the relative 
phase of the condensates. To this aim we focus our attention on the evolution of the operator Q re i, since for the first 
order correlation function in the Rotating Wave Approximation: 



< > oc < exp 



. V^Qrel(t) 



> 



oc exp 



(54) 



for a Gaussian quantum stochastic process (see |26(]), Knowing that the fluctuations of the numbers of the àtoms 
induce phase diffusion (see JÏ3| , Q ) , we can study the influence of the respective fluctuations of the total number 
of àtoms and of the relative number of àtoms on the evolution of the relative phase. In fact the evolution of Q re i 
doesn't depend of P tot as it can be seen in the Eqns. (34a, 34b, 34c, 34d) so that the fluctuations of the total number 



of àtoms can induce diffusion only for the total phase. We will show that the fluctuations of the relative number 
of àtoms will, in counterpart, induce a diffusion of the relative phase. We consider the case when SV is not equal 
to zero, but small enough to make a pertubative development of all the interesting physical quantitics. We assume 



that the traps have different freqüències so that SV = m(u 



,)r 2 /4 



TO&J 2 r 2 /4. 



In fact the real dimensionlcss 



parameter of the perturbative development appears naturally in the calculation and is mii ^^ 2 that we will denote 
v in the following and which represents the ratio between the two characteristic energies of the system. We solve the 



equations (34c 34d) in the first order in v. To find the first order correction for the relative phase operator, Q^, 
we have then to calculate some of the integrals appearing in in its evolution equation with the wave functions given 
at zero order in the hydrodynamic limit (namely ipo+, "00-: 4>o+i </>o-)j some others using the first order correction of 
these functions (in the hydrodynamic limit) and replace the operators by their contribution of zero order in SV. 
In this manner we obtain: 



dt 



1 



(o) (o) 
7+7- 



Nü^ 



(o) (o) 
7+71 



1 



(o) (0) 
7+7- 



P 



(0) 
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^im^fi^V^rM^I^) 



(o) 
7+ 



i(o) 



(55) 



Before calculating the coefficients of the differential equation ( |55[ ) we remark that their scalings with N are different. 
Indeed one of the terms in front of Pj^J is proportionnal to the quantity of the type {(f^?\((>_ ) which is divergent in the 
hydrodynamic limit, and in order to calculate it, corrections to the hydrodynamic limit have to be taken into account. 
The correction of the divergence of the integral is achieved by introducing the healing lenght £ of the condensate 
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defined as, £ = y m $. r ° ) (sec p7||). At the first order in the difference of the traps the sealing of this term with N 
is like vhiN although the others scale as v. Since we are working in the hydrodynamic approximation which is vàlid 
for the great vàlues of N (see p4j), we can neglect the contribution of the terms proportionnal to Qrel m comparison 
with those proportionnal to P^J and in fact forget all the terms which do not have the form of a "00" product. The 
final result is the following: 



QÍ>) 



—v lní— 1 — ; 

5 r J y/N^Nl 



sin(2Aí 



(56) 



where A8q = — B^)/h is the difference between the phase of the zero order wave functions (see (35a , 35b| )) . In 
the first order in v the fluctuations of the relative phase are an oscillating function of frequeney 2À, provided ro is 
constant in time (the latter assumption means that at the time t = the densities of the A and B àtoms have the 
stationary form, and the dynamics is triggered only by the phase difference). It can be noticed that v oc tq and that 
in the hydrodynamic limit we have tq oc N 2 ^ 5 so that v oc A 2 / 5 too. It means that our pertubative developement in 
the parameter v is vàlid only for vàlues of 5V, X and N fullfilling the condition t)<l. At the first order in v the first 
order correlation function is then a ? periòdic function. The physical interpretation of this result is that the relative 
phase between the two condensates exhibits some fluctuations periòdic in time. This dynamical behavior could be 
compared to that of a squeezed state of an harmònic oscillator undergoing free evolution. 

In order to find a diffusion term we have thus to go to the second order in the development. We find (retaining 
only the coefficients presenting a sealing with N as v 2 ln N, where the logarithmic dependance with N can only come 
from integrals involving a product between two functions 0) '■ 



dt 



= -16 



Pi 1 , 5 - 16 



Retó_ 1} | 



,(o) 



) + Re(^° ) |0 (1) ) 



P 



(0) 



Im(^ O) |0- O) )-M 



(57) 



The evaluation of the first order correction of the functions is made in the Appendix. For the final result we have 
only retained the part of the operators proportionnal to Pj^J (assuming that the initial fluctuations with the relative 
phase are small compared to those of the relative number of àtoms) and the part giving effectively the diffusion term. 
We then have: 
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,2^ 



NAN 



nxp?2t + 



(58) 



where A7V = N + — N-. This result shows that at the second order Q re i is growing linearly in time. The first order 
correlation function will then present the usual Gaussian decay in time : 



< W + (t)y-(t) >oc exp 



t 2 



'D 



(59) 



with 



' 25 r 'N + N- V rel 



(60) 



assuming that the mcan valuc of the relative number of àtoms is fixed which means that < P r ( ° } >= 0. 



VII. CONCLUSIONS 

Generelazing the method previously developed in Ref. |Ï4|] , we have made an analysis of the dynamics of two coupled 
condensates adopting a complete field theoretical approach. We have derived the evolution equations of the operators 
associated with the fluctuations of the total and relative number of àtoms and phase of the condensates. The first 
result is that similar ly to the static casc (sec |ï^]), the dynamics of the relative operators is decoupled to that of 
the total operators. Using these evolution equations it is possible to express the equal time first order correlation 
function and show that its behavior is drastically influenced by the difference between the two trap potentials. In 
the case of harmònic traps with equal freqüències the result is that there is no dephasing between the two phases 
of the two condensates. The relative phase fluctuations remain constant in time. But if the traps are different the 
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phases present a dephasing with time. The dephasing process is due to two effects: the first is the non-linearity of 
the dynamics of the system, the second is the usual spreading of a wave packet in quantum mechanics during its 
propagation. We have shown that when the relative number of àtoms is allowed to fluctuate around a fixed mean 
value, these fluctuations induce a diffusion of the relative phase. We have analysed the dephasing time for the case of 

c 2 

two harmònic traps having different freqüències in the limit where the ratio v = m ^ A r ° is small compared to 1 using 
a perturbative approach, and shown that this dephasing time is of the second order in v. It can be expressed as: 

It is worth noticing that the characteristic dephasing time depcnds on the healing lcnght of the condensates, which 
indicates that the phenomenon takes its source in surface effects. 

To conclude we remark that the method developed here could as well be applied to the case of two harmònic traps 
shifted with equal frequency (see pq ) and could be generalized to the case of the mixing of three condensates (see 
for instance and references therein). Into the latter case, the method allows to calculate the spin mixing time of 
different spinor components of the Bose-Einstein condensate. 

We thank Y. Castin and R. Dum for enlightening discussions. 



A. Appendix 

In order to evaluate coeficients in the differential evolution equation of the operator Q re i we have to calculate the 
first order correction to the ip± functions and of the <j>± functions. To this aim we decompose the functions as: 



^+ (r, í) = VA° } 0, t) + A+ (r )e 2lAt + B+ (r )e 



ip_ (r, t) = ?/A 0) (r, í) + A_ {f)e ZlAt + B_ (f )e 



2/Aí 



0+(f,í) = <j>f{r,t) + C + (r)e^ + D + (r)e 



2i\t 



<f>- (r, t) = <£ 0) (r, í) + C_ (r )e ZlAt + D_ (r )e 



2i\t 



-2iAí 



,-2iAí 



,-2iAí 



-2i\t 



(62) 
(63) 
(64) 
(65) 



In this way we take into account only the most relevant Fourier's components of the function. We then introduce 
these forms into the Eqs.(10a 10b) and linearize them in terms of the corrections to the zero order functions. As 
the inhomogeneous terms in the equations for the complex functions A± , B± , C± , D± are proportionnal to SV these 
functions are the first order corrections that we are searching. For the first order corrections of the tp function we 
have the following system of time independant differential equations: 



M 



A+ 
A* + 
B+ 
B% 
A_ 
A*_ 
B_ 
B* 



-SV 



(0) 



(0), 



(0) 



1>. 
^* 






(66) 



where M is the following 8 by 8 matrix: 



M 



AL, 



M_ 



(67) 



with: 



M ++ = 



£j + 2HX 



A^r 2 



C+ + 2HX Nu ^ ] * 2 







Ct - 2h\ 



(68) 
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M__ 



Cq + 2hX 






Cq + 2fi\ Nuo^* 2 



Nu ip 



(0)*2 







(0)2 







Nu ip 





(0)2 



2h\ 



Mi = 













,(°) 















(69) 
(70) 
(71) 
(72) 



where we denoted by £ ,_ the operator £o + uoP+ - with £o = £ + u a p ( ' . In fact as we are working in the 
hydrodynamic approximation we can neglect the result of the action of the Laplacian in C on the first order corrections. 
That means that the above differcntial system becomes in fact an algebraic system of equations. We can write the 
matrix M in the following form: 



M = 

J 
J 



A + M , 
with 



(73) 



where A is the following diagonal matrix: A 



J = 



Proceeding in this way, the formal solution of the system of equations 

1 



-2h\ 








o - 





2h\ 














-2h\ 





. 








-2hX. 



(74) 



|is: 



$(i) 



A + Mq 



8V 



(75) 



where we use a vector notation. In the vinicity of ro all the coefhcients of the matrix Mo go to zero, so that the first 
order correctiontof the ip function have asymptotically the same behavior as the zero order parts. Namely, we can 
write for r — ► ro: 

(76) 

Using the same approach we find the following system for the first order correction to the <fi functions: 










$(o)' 









M 



\c+- 




r Ti i 


c; 




T* 


D+ 






D* + 




r 2 * 


C- 




tí+-Ti 


Cl 




n + _T x * 


D- 




íi+-t 2 


_D*_ _ 




_n+_T 2 *J 



sv 



j,(0)* 
j,(0)* 




where II_| is the operator intcrchanging the signs + and — , and the terms T\ and T2 are: 



(77) 



Ti = 



T 2 = 



NÜ&A+ - 2Nu A + (^^ ) * + 4 0) >f ) - Nu A + (d>^^ + 0<_°>V<_ O) )] - Nu A^^> 

27v mo b;^°VV 0) + nu b*_ (^L 0) + 

2Nu A* + ^U i + ) + Nu A*_(ipf<t>M + 



6<°>>f ) 



,(0) ,(0)* . ,(0)* ,(0)n 

í>vv>+ +0+ ¥>+ ) 



(78) 

(0)n 
+ i 

(79) 
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the integrals of the form Re(^>^ |0=c^) i n the hydrodynamic limit. The formal solution for the first order correction 



We are mainly interested by the first order correction to the <f> functions which will make lead to a divergence in 
f the form Re((/>±^ \4^) in the hydrodynamic limit 
of the system (|77| ) is: 

'[T] - SV 



(i) 



A + M 



(80) 



so that 



where the notations are the same as above. The first vector involves only products of the form ip± tp± 
in the integrand it will contribute in the form ip^ip^ 4*± Vi** whose behavior in the vinicity of the boundary is 



continuous, since the first order corrections to the ip functions behave as the zero ordercounterpart (see ([Tq)). 
r —>■ ro we have: 



$(1)" 













For 



(81) 



which shows that the first order corrections of the <f> function have the same kind of divergence than the zero order 
parts at the boundary, and that the integrals of the kind Re(0£'°' ) Iç^' 1 -*) will scale as v 2 hiN. In the large N limit, 

" (2) 

only these logarithmic divergent terms must effectively be used in the evolution equation of Q rp , ■ Moreover these 
most relevant parts of the first order correction of the <f> functions can be calculated from Eqn. (J8l|) . It then appears 



that the integrals Re( 
instance we have: 



<- (0) ) and Re^ 1 



d 1 : 



) are constant in time, and take part in the diffusion process. For 



(82) 



Careful calculation of these terms allows to derive Eqs.(|58|). 
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